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Abstract 

We study twist-three effects in deeply virtual Compton scattering on an unpolarized spin-1/2 
target. A careful definition of observables as Fourier moments w.r.t. the azimuthal angle allows 
for a clear separation of twist-two and -three effects. Although the latter are power suppressed, 
they give leading contributions to the twist-three asymmetries and do not affect the twist-two 
observables. 
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1. The understanding of the hadron substructure as well as the dynamics of hadron con- 
stituents requires mutual efforts from both the theoretical and experimental side. Theoretically, 
we have to deduce as much information as possible from perturbative and non-perturbative QCD 
calculations. Until now, on the experimental side the nucleon structure has been mostly probed 
in inclusive processes, e.g. measurements of deep-inelastic structure functions in leptoproduction 
experiments, which give access to parton distribution functions. Additional information has been 
gained from exclusive measurements of electroweak form factors sensitive at high momentum 
transfer to the lowest Fock component of the hadron wave function. Since recently it became clear 
that far richer information can be extracted from generalized parton distributions (GPDs) 0, |^, 
which appear e.g. in the processes of deeply virtual Compton scattering (DVCS) ^N f'yN' and 
exclusive leptoproduction of mesons iN — > i'MN'. The GPDs include as limits the conventional 
hadron characteristics alluded to before |[|, 0, 0]. 

The first experimental results 0, ^, ^ show the experimental feasibility of DVCS measure- 
ments. The DVCS process contributes together with the contaminating Bethe-Heitler (BH) one 
to exclusive electroproduction of the real photon and requires detailed theoretical studies. Obvi- 
ously, the first issue is to define appropriate observables and to determine the phase space regions 
that would allow for a clean interpretation of experimental measurements. As has been shown in 
previous studies, diverse asymmetries give an access to the interference term. Here the DVCS sig- 
nal gets augmented by the Bethe-Heitler process, which provides the required handle on the GPDs 
from leptoproduction reactions p|, |^, §]. Certainly, at very large values of the photon virtuality 
the DVCS process is dominated by its leading twist-two approximation arising from the so-called 
handbag diagram (including radiative corrections). Practically the onset of scaling can only be 
judged experimentally as present day theory does not provide necessary non-perturbative tools 
for the computation of soft contributions to such a reaction. For deep-inelastic scattering (DIS) 
it turned out, that the twist-two approximation is valid down to a rather low photon virtuality of 
order of a few GeV. Although the DIS process is given by the absorptive part of forward Compton 
scattering, calculated by means of the same operator product expansion, the DVCS case remains 
an issue for studies. Since the average for all present facilities is about 4 GeV^ or even below, 
it is vital for a clean isolation of the leading twist GPDs, to have a control over power suppressed 
effects. 

Recently, the issues of twist-three effects |10[^ HH and target mass corrections ||20| have been 
addressed in the literature (in the latter case only partially). In the present paper we make a 
step towards accomplishing the goal of handling the higher twist contributions by defining and 
computing twist-three observables in the DVCS cross section with polarized lepton beam and 
unpolarized nucleon target. 
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2. The DVCS hadronic tensor is given by the time ordered product of the electromagnetic 
currents = ip'j^ip sandwiched between hadronic states with different momenta. In leading order 



of perturbation theory it reads |11 



T^.(g,P, A) = 1 1 dxe''-^P2\TUx/2)U-x/2)\P,) (1) 

= -T> a -p llYl. A_ CP p -p i-p p p ] Y2± p ie V — 

F ■ q F ■ q F ■ q 

where we have kept all contributions to twist-three accuracy. Here we have used conventional 

conditions on invariants in the Bjorken limit, — ~ P ■ g = large, = small, ^ = ~q^ /q ■ P = 

fixed, constructed from the vectors P = Pi + P25 ^ = P2 — Pi, and q = (gi + g2)/2. The vectors Pi 

(gi) and P2 (g2) refer here to the incoming and outgoing proton (photon) momentum, respectively. 

The reality of the final state photon implies the presence of only one scaling variable ^. Current 

conservation in the tensor decomposition ([1|) is ensured here by means of the projection operator 

P.. = 9,. - (2) 
qi ■ q2 



This is consistent with an explicit calculation of the amplitude (|T]) to twist-three accuracy |10, 11 
In the above equation the vector V2p is expressed in terms of the vector Vip and the axial- vector 
Alp form factors 

V2p = ^V,p-l^q-V,+ '-'f^A^. (3) 
2 F ■ q 2 F ■ q 

The latter are given as convolutions w.r.t. the momentum fraction x, ® = j\ dx, 

V,p = x) ® v^pix, A2), A,p = x) ® a,p{x, ^, A^) (4) 

of the leading order coefficient functions 

e x) = Q',{1 - x/e - te)-' T{x-^ -x), (5) 

with — (+) standing for the parity even (odd) case, Qi being the electric charge of a quark of 
type-z, and the Fourier transforms of light-ray operators (n^ = 0) being defined by 

' ^x,e,A2)= /^e^-("-^)(P2|^(-«:n)| \ [-k, ^]^{^n)\P^) . (6) 



The general decomposition of the vector and axial- vector amplitudes, in a complete basis of form 
factors to twist-three accuracy, reads 



q-P q-F ^q-F ^q-F ^q-F ^q-F 



q-F q-F '^q-F ^q-F ^q-F ^q-F 
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where = Ap + ^Pp, A^ = iepAPq/P ■ q and the Dirac bihnears are defined conventionally by 
hp = U{P2, S2hpU{P,, S,), ep = U{P2, S2)icjp„^U{P^, ^i), 

hp = U{P2, S2hp7,U{Pu ep = ^2)75^/(^1, ^1), 

with U being the nucleon bispinor. 

In order to simplify our presentation, we introduce a unified convention for the twist-two 
JF = {H, £, H, S} and twist-three = {TC^, . . . , S^} Compton form factors. Replacing the sets 
of Compton form factors JF and JF|. by the sets of GPDs Fi = {Hi, . . . , Ei} and so on, we write 
in analogy to Eqs. (0,^) the decomposition for the matrix elements vip and aip in terms of GPDs 
for quark species i. Here [Hi) and Ei (Ei) are leading twist-two spin non-fiip and spin fiip 
GPDs, respectively, in the parity even (odd) sector. The remaining two sets, containing together 
eight independent functions, belong to the twist-three sector. As mentioned before, the three 
sets of Compton form factors are given as a convolution of the coefficient functions with GPDs, 
generically written as (replace — > JF|. and F ^ F| to get the result in the twist-three sector) : 

^(e,A2)= C^He,x)®F,(x,e,A2). (9) 



In order to deduce the GPDs from the result in Ref. |rT|, we decompose the (axial-) vector 
Dirac bilinears hp (hp) in its twist-two and -three components by means of the Dirac equation 

_ q-h AM^ j ./ q-h q-e\ ~^f A^q-h g-e\] 

' ''g-P (l-e)(A'-ALn)l " V4M2g-P q ■ P J ^ ^ [iM^ q ■ P q-Pjj' 

'q-p^ {i-e){^'-^iA'^\\^M^ )<i-p eq-p) ^ ^ 

+ A- (^!LJL_!L1\\ 

P \AM^q-P q-P) 

where twist-four terms, proportional to gp, have been neglected. The GPDs can now be easily read 
off from Ref. [|r^ . All twist-three GPDs are decomposed in the so-called Wandzura-Wilczeck terms 



(also calculated in Ref. [|T^) and a function F]^'^ that contains new dynamical information 
arising from antiquark-gluon-quark correlations: 

Fl = Ff ^ + Ff. (11) 



According to the analysis of Ref. [ 1 1 1 , the WW parts have the following form 



1 - 1 4/\/f2 



1 4A/f2 
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where the convolution with the PF-kernels 

is defined as previously in Eq. (|) 

W0dF = £^dyw(^^,^Y{y,OF{y,0, (14) 

with the differential operator d{y,C,) = y-^ ~ 

Compared to the results for a (pseudo) scalar target (cf. with Eq. (24) in Ref. |jl2|]), in the WW- 
sector of Eq. (|12|) there appear in addition the functions F^, which arise from the decomposition of 
hp = hp — Ppq-h/q-P and hp = hp — Ppq-h/q-P, respectively. Their calculation is straightforward 
and they read 



4M2 

Ei = ± {^W± ®{H + E)-W^®h], 
Hi = ±U {^-^^^^®H+^W^®{H + E) 



Ei 



4M2 J 

±^{W±®H-^W^(^{H + E)}. (15) 
The antiquark-gluon-quark contributions, 

^ - /; f £ ..i^ [w (-1 -I) -0 . w (I I) |.,(. -oj , 

can be read off from the parametrization of the corresponding operators and result in eight inde- 
pendent functions 

P Pq.p Pq.p ^ Pq.p H Pg.p E 



(see Ref. [[TTl for details). 

It turn out that only the difference F^ — F'^ will enter in the DVCS amplitude, thus, only 
four new GPDs remain at twist-three level. In the WW approximation, i.e. neglecting the quark- 
gluon-quark correlation, all the twist-three Compton form factors are entirely determined by the 
four twist-two GPDs H, E, H and E. 

3. In this section we define the twist-three observables, which are accessible by measuring the 
four-fold cross section for the process e{k)h{Pi) —>■ e{k')h{P2)'~f{q2)'- 



da a^XBy f, 4M2x|\ 



1 



dxBdyd\A'^\d(j) SttQ^ I Q 
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(16) 
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This cross section depends on the Bjorken variable Xb = —q'f/{2Pi ■ qi) (with qi = k — k' and 
qf = — Q^), the momentum transfer square = (P2 ~ -Pi) ^5 the lepton energy fraction y = 
Pi ■ qi/Pi ■ k, and the azimuthal angle between lepton and hadron scattering planes. In the 
twist-three approximation ^ can be replaced by xb via C, = a;B/(2 — xb). Furthermore, we choose 
a frame rotated w.r.t. the laboratory one in such a way that the virtual photon has no transverse 
components^. We fix our kinematics by choosing a negative z component for the virtual photon 
momentum and a positive x component of the incoming electron: k = {E, E sin 6^,0, E cos 6e), 
qi = (g?,0,-|g?|), Pi = (M, 0,0,0) and P2 = (^2, IP2I cos0sin^H, IP2I sin0sin^H, IP2I cos^h), 
where is the azimuthal angle between the lepton and hadron scattering planes. 

The amplitude T is the sum of the virtual Compton scattering (VCS) Tycs and the BH 
amplitude Tbh- The latter one is real and is parametrized in terms of electromagnetic form factors, 
which we assume to be known from other measurements. The azimuthal angle dependence of each 
of the three terms in = |Tb//P + iTycsp + X, where X = Tvcs'^bh + '^vcs'^bh, is given in our 
frame by a finite Fourier sum. In the case of an unplolarized or longitudinally polarized lepton 
beam, the interference term X and the squared DVCS amplitude iXv/c^p may be written as 

f A2 2 



XBZ/^PlP2(-A2) O^Z^^ 



cos(m</)) + Xs^ sin(m(/ 



ITwcs P = ^ {C^^ + K [cr^s cos(0) + A .r^s sin(0)] } . (17) 

The coefficients cf, as well and sjj'^*"^ arise at the twist-two level and their dependence 

on GPDs has been calculated in Refs. [0, |^. cj, cf , sf, cf^^^, and sf^*"^ provide an additional 
angular dependence and are given in terms of twist-two and twist-three GPDs, while the termsQ 
that are discarded here are either proportional to cos (30) [cos (20)] or sin (30) [sin (20)] for the 
interference [squared DVCS] term. Note that all c's and s's are independent. There is an 
important difference between the interference term and the squared DVCS amplitude. The former 
has an additional dependence due to the lepton propagators^: 



Q^Vi = {k- q^f = Q2 + 2/e ■ A, Q^V2 = {k - A)^ = -2A: ■ A + A^ 



with 



■ A = -^h + 2ir COS0 - — (1 - xb(2 - y)) - -^(2 - y)xl\ {l + O {m'/Q\ AVQ')} • 



^It is a reference systeni related to the centre-of-mass system of Q by a boost of the hadron in the z-direction. 
^They are induced at twist-two level by the gluon transversity, which is perturbatively suppressed by ag and is 
contaminated by twist-four contributions. 

•^For convenience we scale them with respect to . 
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The 1/Q-power suppressed kinematical factor, appearing also in the Fourier series ([T7|), 

^=|-|^(l-^B)(l-y)(l-^)} ^ {l + 0(MVQ^AVQ2)} (19) 
vanishes at the kinematical boundary = A^^j^^, determined by the minimal value 

Since X and \Tbh\'^ + |^c5p are charge odd and even, respectively, those experimental fa- 
cilities that possess beams of both charges can separately extract the BH-DVCS interference 
term and the sum of the squared BH and DVCS terms. Combining the charge asymmetry with 
different nucleon/lepton polarizations and azimuthal asymmetries p, |^, || one can extract, in 
principle, the real and imaginary part of the four twist-two amplitudes 7i, Ti, and S from 
the interference term Since I is linear in amplitudes involving GPDs, one gets access to the 
twist-two GPDs, which, however, are convoluted with the real or imaginary part of the coefficient 
functions. The precise definition of the Fourier coefficients is crucial for the interpretation of ex- 
perimental results, since the w-channel lepton propagator gives a contribution that behaves like 
Vi = [1 — y — K coscj) + (9(1/ Q^)] /y. For instance, in the case of the single lepton spin asymme- 
try it may be possible that in some kinematical regions, i.e. large y and not too small A^/ Q^, this 
propagator effectively induces a sin (20) term that is canceled by a sf sin (20) term resulting in a 
(fake) sin dependence of the cross section. 

The procedure, outlined above, can also be used to access the twist-three contributions. For a 
successful separation of definite twist components, one should compensate for the strong azimuthal 
dependence of the lepton propagators in the BH amplitude. The cos(m0), sin(m0) components 
are e.g. unraveled by weighting the cross section with V1V2 {cos{m(f)) , sin(m0)}. If this is not 
done, the magnitude of the twist-three effects can be judged by the distortion of the leading twist 
angular dependence {cos(0), sin(0)} ^Pj^^. 

In the following we only present the Fourier coefficients to twist-three accuracy for an un- 
polarized nucleon target. A complete analysis will be given elsewhere A straightforward 

computation provides the following analytical results for the squared DVCS amplitude 

CSf = 2(2 -2y + y-)CZ''' (-^, ^ , (20) 

DVCS ^ g 2 — y p ^DVCS f -pcS 
,DVCS _ Q y A-DVCS f TcS 

2-xb' 



Both twist-three coefficients cf^R^ and are given as interference of twist-two GPDs with an 
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'effective' twist-three GPD 

^^fif = -2^ [yt^-^ + " •^-) ' ^^^^ 

where are defined in Eqs. ( p!2[jT6| ). Surprisingly they have the same functional dependence as 
the leading twist-two function 



Cp^^ (-^. ^*) = TTTz — ^1 ~ ^b) [m* + nw) - xl ( HE* + £n* + ns* + ^ ti:* 

(z — XbJ I 



'^B + (2--B)^j§Jj)f£--x^^fr}. (22) 



For the interference term we found the same property for the cos(20) / sin(20) coefficients: 

4,unp = -8(2 - y) |(2 - xb)(1 - y) - (1 - xb)(2 - yf [l - ^) | Re^^ (^) 

+ 8(2 -y){l- y)xB{F, + F,)Re\-^{n + S)+n], 

1 2 - xb J 

cjunp = -8(2 -2y + 2/2)ReCj„p (^) , sj,,p = 8y(2 - y)lmCl^ (^) , 

cjunp = -le^ff Re Cj„p (^^^^) , = 16-^Im C^^ {j^^^) . (23) 

Here the twist-two function 

Cp (^) = F,n + ^^(^1 + F,)n- ^^F,£ (24) 

also depends on the Dirac and Pauli form factors Fi and F2, respectively. 

Note that twist-three GPDs, having generical discontinuities at = ^, enter in the Compton 
form factors in a singularity free combination just as for a (pseudo) scalar target. The convolutions^ 
of the coefficient functions Cf with the IV-kernels make this property transparent: 

^•«) = TT^^ + {|^-'^'«-^' ^ + ^, _ _ ^^f tia (25) 

where we have used a new convention for the coefficient function 

Cf^\i,x) = 7^ In T {x - -x}, (26) 



*This can be done by means of J^^ ^ ^-x-io ^ (|' f ) = ^f^^ ( c-y^io ) ^^^^ proper symmetrization. 
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with the — (+) sign standing for T = {H(H), S{£)} and 



4M2 
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We should note that the kinematical factor in front of JF-*- drops out in the final results p0|j23| ) 
and, therefore, the 1 — A^j^/A^ dependence in the Fourier series ( jTTD will not be altered. More 
precisely, the dependence in the 'effective' twist-three amplitude (^) generates the following 
terms: 

C?np [^''^) = C?np -^^ = 0) + -^(F, + F,) E Cf ^^ ® H^i (28) 
(^^^ = C^p^^ (^^^ -F*; = 0) 

+ 2xB (^E ® (^^ + Hi)^ ((2 - ^b)7?* - xb{s + ny) 
+ 2xB ( E cf^^ ® ^0 ((2 - ^b)(h + ^)* - XB{n + sy 



4. In this paper we have defined twist-two and twist-three observables in such a way that 
the twist-three contributions, suppressed by (9(A_|_/Q), do not alter the leading twist angular 
dependence. Thus, the corrections to the twist-two Fourier components are expected at twist- 
four level (9 (M^/Q^, A^/Q^). Combining charge and spin asymmetries together with a careful 
extraction of the different Fourier coefficients allows a separate measurement of all coefficients of 
the interference and the squared DVCS term. In the latter case, the subtraction of the squared BH 
amplitude from the data is necessary. At the same time it drops out in single spin asymmetries. 

For a lepton beam with definite charge, a separate measurement of the interference term is 
not possible. In this case the leading twist cos(0)/ sin(0) dependence in the cross section, which 
stems from the interference term, gets corrected by the twist-three contribution of the squared 
DVCS amplitude. Although such corrections are kinematically suppressed by xb^^A^/Q^, their 
numerical value depends on the magnitude of the form factors JF. 

Let us also mention that for a (pseudo) scalar target the Fourier coefficients fulfill certain 
constraints derived from general arguments [T^. They provide a basic test of the applicability of 
the operator product expansion at a given Q-scale. We expect similar constraints for a spin-1/2 
target, however, their derivation will require an evaluation of the squared amplitudes for longitu- 
dinally and transversely polarized targets. Of course, one can not test the dominance of leading 
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contributions in this way or the magnitude of higher twist and as suppressed contaminations. 
This can be judged by comparing the measured scale dependence of the coefficients q and Sj with 
the ones arising from the evolution of GPDs and power suppressed contributions in a sufficiently 
large Q^-interval. Note that this task is complicated for experiments with fixed center-of-mass en- 
ergy, since the kinematical variables xb, y and Q do not form an independent basis in this setting. 
On the theoretical side the complete perturbative NLO corrections have been calculated while an 
estimate of twist-four effects still has to be done. 

Extraction of the leading twist-two and -three Fourier coefficients allows to test models for 
GPDs and even to measure them in single spin asymmetries on the diagonal^ x = ^. Having tested 
those predictions and processing to higher twist, one can study the validity of the 'Wandzura- 



Wilczek model' for twist-three functions suggested in Ref. [TT|. A significant deviation from the 
latter would imply an essential contribution from antiquark-gluon-quark correlators. 

To conclude, the full experimental exploration of the deep electroproduction of real photons 
with leptons of both charges and polarizations will lead to a direct confrontation of the data with 
theoretical predictions and will result in systematic tests of our understanding of the quark-gluon 
content of the nucleon via generalized parton distributions. 
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